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Abstract 

A method is described for levitating micron-sized few layer graphene flakes in a quadrupole ion trap. 
Starting from a liquid suspension containing graphene, charged flakes are injected into the trap using the 
electrospray ionization technique and are probed optically. At micro-torr pressures, torques from circularly 
polarized light cause the levitated particles to rotate at frequencies >lMHz, which can be inferred from 
modulation of light scattering off the rotating flake when an electric field resonant with the rotation rate 
is applied. Possible applications of these techniques will be presented, both to fundamental measurements 
of the mechanical and electronic properties of graphene and to new approaches to graphene crystal growth, 
modification and manipulation. 



1 Introduction 

Since its discovery in 2004 [1 , graphene has received a tremendous amount of attention directed both towards 
understanding its fundamental properties and seeking applications for this new material [2][3]|1]. From the 
physics perspective, graphene is the first truly two dimensional (2D) system, with electronic, mechanical, and 
thermodynamic properties all determined by the structure of a single sheet of carbon atoms. It seems likely that 
the revolutionary applications for graphene in the future will make use in some way of the unique properties of 
an intrinsically 2D system. 

While graphene is an ideal 2D material, it nonetheless must be coupled to the outside world in real experi- 
ments and, in practise, attached to a 3D substrate in some way. By now it is known that the substrate can limit 
mobility of graphene electrons [5|[6]fr, and consequently the recent experiments demonstrating the fractional 
quantum Hall effect in graphene ^ ^ were conducted on samples that were suspended (i.e. the substrate was 
etched away) in the region of measurement. Locally suspended samples were also used for the first measurements 
of the mechanical properties of graphene [TU] [TT] , but it is probable that loss mechanisms and residual stresses 
in these experiments are determined by the substrate |12| . The thermodynamic measurements of graphene are 
perhaps most likely to be hindered by the presence of a substrate, since the expected melting temperature of 
graphene is in excess of 3000 K, higher than that of any other material. 

It is possible to avoid coupling to a substrate altogether if graphene is levitated using the particle trapping 
technologies that have been perfected in recent decades. Indeed, optical trapping of graphene suspended in solu- 
tion has recently been demonstrated [T3] . Diamagnetic trapping is possible for graphite [I^ . For measurements 
of graphene, however, quadrupole (ion) trapping has several advantages: it is compatible with ultrahigh vac- 
uum (UHV), low, and high temperature measurements; it can provide tight confinement of very small graphene 
particles; finally, ion trapping has the advantage that charged graphene flakes will tend to remain flat when lev- 
itated due to electrostatic repulsion. Interestingly, ion trap techniques have been applied previously to graphite 
particles for studies of interstellar dust [H] [S] [IT] [TH] ■ Ion trapping of micron-sized particles has also been 
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developed as a test-bed for quantum information processing techniques [HJIIO] and for materials studies of 
small particles HJE^. 

Below, I describe a quadrupole trap optimized for graphene measurements and present preliminary data. 
Graphene flakes originally suspended in liquid are injected into the trap vacuum chamber using the electrospray 
ionization technique [231. Trapped flakes are detected optically. The optical absorption and low mass of the 
flakes mean that the particles absorb angular momentum from circularly polarized light ^24) (25. and begin to spin 
rapidly at low pressure (p =l/itorr). This spinning is directly demonstrated by the observation of rotational 
resonance - the modulation of optical scattering from the particle at a well defined frequency of an applied 
electric field - at frequencies above 1 MHz. This high rotation frequency, facilitated by the ability of graphene 
to withstand centrifugal tension during rotation, is, to the author's knowledge, the largest ever measured for a 
macroscopic trapped object [24 ] [26] p 7 ] [27] . 

The experiment remains to be optimized, and all the data presented below are likely taken on multilayer 
flakes, rather than on single layer graphene. Nonetheless, it is hoped that further refinements will enable 
levitation of single layer graphene and that rotation resonance measurements will effectively capitalize on the 
powerful techniques of analysis developed for magnetic resonance and ultimately provide a wealth of information 
on the levitated sample and its internal properties. 

I provide below a detailed description of the experimental design, including the development of many formulas 
necessary to understand the requirements for trapping graphene, the expected optical signal from the trapped 
sample, and the torques that will affect graphene motion in the trap. After a presentation of the data, I 
conclude with a consideration of possible improvements to the current design and a discussion of applications 
for graphene trapping, both for fundamental measurements and as a possible new environment for graphene 
crystal modification and growth. 

2 Description of the Quadrupole Trap 

At its simplest, a quadrupole trap is comprised of two or more electrodes in a configuration where the electric 
field £^=0 at some point r = in space away from the electrodes [lH][in|- If only DC potentials are applied to 
the electrodes then confinement of a charged particle at the point r = is always unstable, but stable trapping 
is possible with AC applied fields under appropriate conditions. The equation of motion is: 



where q is the particle charge, m is its mass, r2t/27r is the frequency of the trapping field, and dots above variables 
are used here and below to signify time derivatives. 7i, is the velocity damping rate, which will be extremely 
important for graphene flakes, due to their small mass (typically '-^ lO"^'^ kg for a /im-sized monolayer) and 
large surface area to mass ratio. 

While it is possible to get exact results valid in the limit r — > by linearizing Eq. 1 [30], a more intuitive 
picture of ion trap dynamics comes from the pseudopotential approximation [25]. In this approximation the 
particle is assumed to be rapidly oscillating and experiences a spatially nearly uniform electric field during each 
oscillation period. The position of the particle is separated into an oscillatory "micromotion" S and a position 
R of the particle averaged over an oscillation period: r = R + S. Trapping occurs because in an oscillating 
nearly uniform electric field in the absence of damping, d and d (oc force) are 180° out of phase: a particle in the 
neighborhood of r = will experience a force towards r = when it is furthest away from r = and away from 
r = when it is closest to it. If £^ — )■ at r = 0, then the force on the particle averaged over its micromotion 
will be directed towards r = 0, leading to stable trapping. The effect of finite damping is to reduce the phase 
difference between S and 5, which will tend to reduce the average trapping force directed towards r = 0. 

The above arguments can be used to derive an equation of motion for the particle position averaged over 
the micromotion |29j[16j: 
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Figure 1: Cross section of the apex of the trap electrodes (grey regions) and the surrounding pseudopotentiaL 
The trap is made from two coaxial pieces of conically tapered stainless steel. The oscillating trap voltage is 
applied to the outer electrode, while the inner electrode is held near ground. The surrounding vacuum enclosure 
(not shown) is also grounded. The tip of the inner electrode extends beyond the face of the outer electrode by 
about 200 /im. 



From the formula for the pseudopotential, it is seen that a smaller Qt leads to greater confinement (provided 
the pseudopotential approximation remains valid), and that the confinement rapidly weakens when > Qf 

The trap used in the experiments described below is similar the stylus trap developed for ions 3l|. It is a 
coaxial arrangement of two conically shaped pieces of stainless steel, the apex of which is shown in Fig. 1. The 
AC trap voltage is applied to the outer electrode, while the inner electrode is kept near ground. Additionally, 
the chamber surrounding the trap (^2 cm away from the apex) is also grounded. Calculated |32j values of 
E^{R) are plotted, showing the pseudopotential minimum point about 0.8 mm away from the trap apex. 

Along the axis of symmetry of the trap in the neighborhood of r = 0, E{r) may be approximated as: 

E,^-^Vout and E.^Ey^O. (4) 

Here, Vout is the amplitude of the voltage applied to the outer electrode, and zq is a parameter determined from 
modeling the electrode configuration, and is 1.84 mm for the trap design depicted in Fig. 1. When '^y = 0, Eqs. 
2-4 are readily solved to get oscillating solutions with: 

_ I q Vout 
''^-V2mn,zl 

Similarly, for radial motion away from the axis of symmetry, oscillatory solutions can be obtained with uj^ = 
LOy = Wz/2. These solutions are obviously only valid in the regime lOz ^ ^t- However, the trap will remain 
stable when w^/^t <0.32[l9]. 

Nonzero biases applied to the inner electrode on the trap will produce an electric field at r = Q of Ez = 
—Vin/zi, where 2;i=6.4mm is another parameter that can be determined by trap modeling |32| . Nonzero l^n 
at DC shifts the trap minimum position along the z axis, while l^n applied at finite frequencies can be used to 
determine lUz (and thus q/m) by observing resonance behavior. 

Quadrupole traps are usually designed to maximize confinement while maintaining stability. For atomic 
systems, this requires that fit > 10^ see" ^, and excitation at this frequency is typically provided with a 
tuned resonator. For graphene particles charged using electrospray ionization, q/m^lO-lOO C-kg~^(10~^ — 
10~^ |e|-amu^^), much smaller than the values for atomic ions. Consequently, the trap excitation is conveniently 
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Figure 2: Predicted velocity relaxation rate, 7„, of graphene as a function of pressure. 

provided by a circuit |33j using a high voltage op-amp j34j . This allows for a maximum voltage amplitude of 
400V at a maximum frequency of about 100 kHz. Using Vout=300V, r2t/27r=30 kHz (values typically used in 
these experiments), and q/ni=10C-kg^^m Eq. 5 for the trap in Fig. 1, I obtain: =3300 see" ^. Thus, the 
assumptions of the pseudopotential approximation are well satisfied during the experiments discussed below. 

When a particle is confined by a quadrupole trap, fluctuations about the pseudopotential minimum will 
ultimately be determined by Brownian motion 13! . Using the equipartition theorem: 

^-kBT^^-mcoliz^). (6) 

For temperature r=300K and to = 10"^^ kg, (z2)i/2= 20 /im for the hypothetical particle discussed above. 
This dimension is small compared to the size of the trap, but it is readily measurable, and can be used to 
provide information about the trapped particle. Finally, the acceleration of gravity, g, will displace the particle 
from the trap center by a distance g/w^, which is a few /im in this experiment and is generally not observable. 

3 Particle Damping 

From very crude kinetic theory in the free molecular (Knudsen) regime, the velocity relaxation time of a particle, 
1/7.1,, is the time it takes to collide with its own mass of the gas molecules in its surroundings. For single layer 
graphene, this time is just the time it would take for a monolayer of gas molecules to attach to the surface 
(if they all stuck and had a mass comparable to C atoms). The pressure-time product necessary to deposit 
a monolayer is known to surface scientists as the Langmuir unit =10~^ torr-sec. Consequently, for a graphene 
monolayer p/jv — 10~^ torr-sec. More refined kinetic theory |35| gives a very similar result for graphene in a 
N2 ambient at 300 K. An important point is that for a 2D object like graphene, the mass and collision cross 
section are proportional to one another, so jy is independent of the lateral size of the graphene layer. However, 
The number of layers, n, will affect the mass, but not cross section. Thus, measurement oi p/jv can be used to 
determine n. 

Using these ideas 7^ versus pressure is plotted in Fig. 2 for graphene for n = 1 and n = 10. As mentioned 
above, finite 7„ weakens the trapping pseudopotential in Eq. 3 unless 71, < ilf For the trap parameter Qt — 
200,000 sec~^, this requirement means that the trap can only operate at pressure below 1 torr. On the other 
hand, finite 7^ is necessary, since - without a drag force - particles could never become bound in the trap. 
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Consequently, the experiments discussed below are conducted by trapping the particles at around 0.5 torr and 
subsequently pumping the system down to 10"'^ — 10~^ torr for most measurements. 

While the picture above is adequate to motivate the design of this experiment, it is worth noting that many of 
the assumptions underlying simple kinetic theory may be unwarranted for a 2D object in the Knudsen regime. 
For spherical objects the values of 7i, determined assuming purely specular and purely diffuse scattering of 
molecules off the surface differ by only about 15% [3S]. However, for a 2D plate moving in a direction parallel to 
its surface, molecules specularly reflected off the surface do not transfer momentum to the plate in the direction 
of its motion, and thus would not contribute to 7^. Purely specular scattering would also not contribute to 
7ai the angular velocity relaxation rate for a plate rotating on an axis perpendicular to the plate. The relative 
contributions of diffuse and specular scattering off graphene is unknown, but it is notable that crystalline 
surfaces with a high degree of perfection are being developed as "atom mirrors" in which the ratio of specular 
to diffuse scattering is significant [36; . 

4 Generation of Charged Flakes by Electrospray Ionization 

Graphene flakes are introduced into the trap using the electrospray ionization technique. Originally developed 
for mass spectrometry |23j . electrospray has also been applied to inject micron-scale charged particles into 
quadrupole ion traps |21j.|19j.|20|. In the electrospray technique, a liquid suspension containing the particles is 
ejected from a capillary tube held at high voltage. Ejected droplets shrink in size due to a combination of liquid 
evaporation and droplet fission, until only dry charged particles remain suspended in the chamber gas. The 
large voltages required for electrospray mean that it is most easily performed at or near atmospheric pressure, 
with charged particles subsequently introduced into a higher vacuum environment through a pinhole orifice 
located near the electrospray emitter tip. 

The suspensions used in these experiments are similar to those developed by Hernandez et a^[37][5S]- These 
workers prepared suspensions by ultrasonication of graphite flakes in a variety of liquids. The resulting mixture 
was then centrifuged to deposit coarse material, while finer particles remained suspended. The particles that 
remain in suspension were shown to be micron-scale flakes of few layer graphene, with n typically in the range 
1-10. 

In order to use these suspensions for electrospray, it is extremely important to minimize non-volatile im- 
purities in the liquid, since these impurities will accumulate on the graphene flake during liquid evaporation, 
and can ultimately outweigh the residual graphene if care is not taken. While many liquids are effective for 
graphene suspensions [37] I have chosen to use an isopropyl alcohol (IPA)-water mixture (3:1 volume ratio), 
since IPA and water are commonly used in mass spectrometry and are available in "MS grade" (low in ionic 
contaminants and <lppm residue after evaporation) ^39^ . Care must also be taken to ensure that only clean 
and inert materials are brought into contact with the suspension and that it is not contaminated with foreign 
particulates. For these reasons I give a rather detailed description below of the materials and techniques used 
to prepare and deliver the suspension to the emitter tip. 

The suspensions are prepared in 8mL glass vials with PTFE-faced phenolic caps (VWR #14230-824). After 
the vials are rinsed in IPA, 5mg of graphite flake (Aldrich #332461) are weighed and introduced into the vial. 

5 mL of liquid is then added, and the cap is closed. The vial is then placed in an ultrasonic bath (Branson #1510) 
in a position so that the fluid level in the bath matches that inside the vial, and the sample is ultrasonicated 
for 30 min. All operations where the suspension is exposed to air are performed in a clean room. 

After ultrasonication the vial is centrifuged (Drucker #642E, approximate acceleration^lOOO^) for 30 min. 
After centrifugation the supernatant has a typical optical attenuation coefHcient (measured at 532 nm) of 80 m~^. 
Using the attenuation dependence on concentration determined by Hernandez [57] of 2500 L-m~^-g^^, the typical 
concentration of graphene flakes in the suspension is about 30mg-L~^. 

The supernatant is drawn through a syringe needle (Hamilton #90122) into a glass syringe (Hamilton #1750, 
0.5 mL volume, with PTFE luer lock connection), which is placed into a syringe pump (Harvard Apparatus, 
Nanomite) located at the top of the quadrupole trap vertical column (Fig. 3). Upon exiting the syringe the 
suspension travels through PEEK 1.6 mm diameter tubing and fittings to the stainless steel electrospray emitter 
(New Objective Corporation, 100 /im tip inner diameter). The metal tip is connected to a wire that allows it 
to be biased at high voltages. 

During operation the emitter tip is surrounded by N2 at atmospheric pressure and is about 6 mm above a 
25 fim diameter grounded stainless steel pinhole aperture (Edmund Optics). Because contamination may arise 
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Figure 3: Diagram showing the particle injection system and the gas handhng of the experimental apparatus. 
The minimum obtainable pressure in the chamber is <l/Utorr 



from particles in the gas surrounding the emitter tip that come into contact with exiting charged droplets, 
the N2 is filtered (Swagelok, SCF Series, 0.003 /xm) prior to entering the chamber surrounding the tip. The 
chamber walls are made of glass, so electrospray emission is visible when the tip is illuminated with a laser 
pointer. Onset of electrospray emission typically occurs at Vtip = ±2kV, with optimal performance around 
3kV. Optimal flow rates during electrospray are typically l-3/iL-min~^. During particle trapping a roughing 
pump evacuates the chamber downstream from the pinhole to about 0.5torr. After a particle is trapped, a 
gate valve between the electrospray source and the main chamber is closed (Fig. 3), allowing the chamber to 
be pumped with a turbo pump to below 10~^ torr. For measurement of pressures > 10~^ torr a capacitance 
diaphragm gauge and a convection gauge located close to the trap center are used. A cold cathode ionization 
gauge is located downstream from the trap to measure lower pressures. Close proximity between the trap center 
and the ionization gauge was avoided to prevent possible discharging of trapped particles, and consequently 
there is considerable uncertainty in the pressure measurements below p < 10~^ torr. Finally, an N2 variable 
leak is located downstream from the ionization gauge to enable control of the chamber pressure. 



5 Optical Apparatus 

During and subsequent to trapping, particles are imaged from light scattered at small angles from the direction 
of an illuminating laser (Fig. 4). The light source is a A =532 nm wavelength laser (Lasermate Corp.) with 
adjustable output power. Laser power is levelled during measurements by using a photodiode and a PID 
controller. Power flux at the center of the beam is determined by scanning the laser across a (non graphitic) 
trapped particle and fitting the scattered intensity to a gaussian. The power flux calculated in this way is 
3300W-m-2.mW-i. 

Care must be taken to avoid light scattering into the collection optics that could wash out the signal. For this 
reason, coated windows (Lesker VPZL-450AR) are used for the optical entrance and exit to the trap chamber. 
Additionally, the trap outer electrode is conically tapered to minimize scatter from corners and edges. Also, 
laser line filter is placed in the collection optics to minimize; intcirfcrcncc from ambient light. 

The collection optics consists of two 40 mm diameter lenses located outside the trap vacuum chamber, 
120 mm away from the trap center, with magnification factor of 4.2. A beam splitter is used so that trapped 
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Figure 4: Diagram of optical instrumentation for the experiments. The electric field of the light originating 
from the laser points out of the diagram. Light received at the detectors on average has scattered 20° from the 
incident beam direction. 



particles may be simultaneously viewed by both a low light level CCD video camera (WATEC 120N) and a 
low noise photodiode (FEMTO FWPR-20-SI) with a 20 Hz bandwidth. A 500 fim diameter pinhole aperture is 
positioned in the focal plane in front of the photodiode to minimize the contribution of background light to the 
detected signal coming from trapped objects. Finally, a 1/4 wave plate on a rotation stage allows for control of 
the light polarization illuminating the trapped particle. 



6 Interaction of graphene with electric fields 

The graphene samples studied in this experiment are irregularly shaped flakes whose size will in general not 
be small in comparison to the 532 nm laser wavelength. However, to get a rough idea of the expected signal 
from scattered light, I will estimate the interactions of a graphene fiake with an applied electric field in the 
quasi-static dipole approximation. In this approximation, the dipole moment, p — eoaE, is estimated from 
the polarization induced by a uniform electric field, a is the polarizability tensor, and eo is the permittivity of 
vacuum. Analytically tractable results for a circular disk may be derived by determining the polarizability of 
an oblate spheroid |40| with semimajor axis, a, and allowing the semi-minor axis, b, to go to zero. Inside the 
spheroid, the complex permittivity is: 



e = eo 



(7) 



where ctsd is the volume conductivity of the material and lo is the angular frequency of the applied field. For 
an oblate spheroid, an , the response to an electric field oriented in the plane of rotational symmetry. En , is: 
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where /? is a geometical factor = 7r&/(4a) for <C a |3D]. To obtain an expression for a disk, I equate 0-20 
with the product of 1T3D and 46/3, the mean thickness of the spheroid in the direction of its axis of symmetry 
averaged over its cross section. The result is: 
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Figure 5: Predicted heating (versus laser power) and cooling (versus time after the heating source is turned off) 
of a graphene single layer in vacuum. Heat is assumed to be removed from the flake by black body radiation into 
a surrounding environment at T = 0, and the heat capacity is assumed to be Sfce per C atom. The laser beam 
diameter is 0.5 mm, and 2.3% absorptivity and emissivity are assumed. Size effects and temperature dependence 
of the heat capacity have been neglected and will become increasingly important at low temperature. 

6.1 Optical Scattering and Absorption 

At optical frequencies it is now well established experimentally [41][42][3] that CT2d — for monolayer 

graphene. Using this value and a; = ojl = 3.54 x 10^^sec~^(for A=532 nm) yields a value of 4 x 10^'^ for the 
term in brackets in Eq. 9 with a— 0.5 fim. Thus, at optical frequencies for micron-scale flakes with one or a few 
layers: 

an = lira . (10) 

This value will scale linearly with n, the number of layers in a flake, provided n is small. Using this approximation 
the optical absorption cross section, aabs, of a flake oriented perpendicular to the incident radiation is [40] : 

27r 

(Jabs — — = mra'^ x = 0.023 x nna^. (11) 

A " Aeohc 

This picture predicts that a graphene monolayer absorbs 2.3% of the light impinging on its area, in agreement 
with experiments |41| . This absorption has significant implications for optical measurements, since sample 
heating can be large in a high vacuum environment at laser powers above a milliwatt (Fig. 5). 

To determine the visibility of a flake I next calculate the forward scattering differential (per solid angle il) 
cross section of a flake oriented perpendicular to the incoming radiation using the same assumptions |40j : 

da\ „, / 1 \ ^ /27r \ ^ , ,o / 0.023 x nna 



[t) ni-^^ - ) • (12) 

The quantity mra^ is proportional to the mass of the flake. The result of Eq. 12 is plotted in Fig. 6 as a 
function of particle mass using typical experimental parameters. These calculations show that a single layer of 
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Figure 6: Predicted optical signal and forward scattering cross section for graphene as a function of its size. 
Finite size corrections to dipole scattering are neglected. 



graphene with a 1 /imxl /xm area should be visible in the experiments. It is important to note, however, that 
visible objects will not conform to the condition that a ^ A. This, combined with the fact that flakes will be 
irregular in shape (and possibly contain regions with different layer thickness), means that these results are at 
best rough estimates of the expected optical scattering from trapped flakes. 



6.2 Optical Anisotropy 

The results above were all derived for a flake in a plane oriented perpendicular to the direction of the incoming 
radiation, when E always lies in the plane of the flake. Because of the thinness of graphene (and the anisotropy 
of graphite in general), ay ^ a± ~ 0. Consequently the optical response of trapped graphene should vary with 
time if the flake orientation changes with respect to the incident radiation. If the flake is randomly oriented in 
linear polarized light (LPL): 

7 \ mean r, / i \ max / , \ min 

da \ 2 / da \ , I da\ , , 

= o 3?^ and — =0. (13) 



^dQ J ^pj^ 3 \dfl J \dQ J 

The maximum value is that determined from Eq. 12, and will be the same for linear or circular polarized 
incident light. Because circular polarized light (CPL) effectively averages over orientation around the direction 
of incident radiation: 

^^^mcan 3 / da \ / da \ _ 1 / da \ '""^ ^^^^ 



= - 



dri/(-,pL 3 \dfl J \dfl J fjpj^ 2 \dQ J 

These equations predict that temporal fluctuations should approach 100% in LPL and 50% in CPL. The slow 
(20 Hz) response of the photodiodc in the experiments means that these fluctuations will be averaged out at 
high pressures, however. 



6.3 Optical Torque 

Because electromagnetic waves carry angular momentum as well as energy, absorption of light (Eq. 11) must 
also convey torque, N = p x E, to the particle if the incoming radiation is circularly polarized: 
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Nj_ = Ij_ujj_ = ± — X 0.023 X n7ra^ (15) 

where S is the incoming radiation power flux. N_i_, wj^, and Ij_ are respectively the torque, angular velocity, 
and moment of inertia associated with motion around an axis perpendicular to the plane of the flake, which 
is assumed to be oriented perpendicular to the direction of incident radiation. For a circular flake with uni- 
form density: I± = ma^/2 = 7rna^p2D/2, where p2B is the 2D mass density of single layer graphene=7.6 
xlO^'^kg-m-^. Thus: 

S 2 

tj_L = ±— X 0.023 X . (16) 

Note that since both torque and mass are linearly dependent on n, n does not appear in this formula. For 5 mW 
laser power and a =0.5 ^m, ui± = 1.1 x lO^sec"^. In circularly polarized light, the angular rotation velocity 
of a graphene flake will increase until the optical torque is matched by frictional drag at Wmax = ^±/ju} This 
result implies that at low pressures, when <lsec~^, rotational frequencies above 1 MHz are possible. The 
thermal rotational angular velocity is given by: 

lfcBT^^/a.^ (17) 

which, for a 1 /im diameter single layer graphene flake at BOOK, leads to = 2 x lO^sec^^. Thus, at low 
pressures, light-induced rotation will exceed thermal rotation of the graphene flakes. 



6.4 Low Frequency Torques 

At frequencies <1 GHz, the bracketed term in Eq. 9 3>1 and: 

16 , 16 



( 1 



(18) 



The real part of this expression is simply the static electric polarizability of a thin disk of radius a |43, . This 
term produces a torque along an axis in the plane of the flake when E is at an angle 9 from the flake plane: 

iV|| = Ej_eoa\\E\\ = — coa^ E'^ sin 9 cos 9 = (19) 
3 

For 2D objects the moments of inertia have the property: I^^i + I\\2 = I-Ly ^^nd for a circular uniform disk: 
/|[ — ma? /A. Thus, for a disk: 

wil = ^—f^^E'^smecosO. (20) 
" 3 m 

Typical low frequency electric fields that can be conveniently applied in the trap (by applying a voltage to the 
inner electrode) are ^ 1000 V-m"-'^. This value is comparable to the electric field that the particle experiences 
due to thermal motion away from the trap minimum. For a=0.5 /im and m = 10^^^ kg, tl;|| is at most 5x lO^sec^^, 
a value that can be considerably larger than the value for the angular acceleration from optical fields calculated 
above. 

Finally, the imaginary term in Eq. 18 leads to a braking effect, since dissipation occurs when charge moves 
on the spinning flake to shield external electric fields oriented in the flake plane. For a circular disk: 

^UE,)^u.J..= '^^4El (21) 

At low frequencies, it is inappropriate to use the optical value for cr2D, since the conductivity will depend on 
the number of carriers in the flake (and possibly on quantum size effects in a mesoscopic system). However, 
to obtain an order of magnitude estimate of the braking rate, I use (T2d = e^/Ah, a =0.5 /im, m = 10^^^ kg, 
and i?=1000 V-ni^^to get 7aj(-E||) = 10~^sec~^. From a perusal of Fig. 2 it is seen that juj{E\\) is neghgible 
compared to damping from background gas (assuming jy = 7^) for the conditions of this experiment. However, 
damping either from applied electric fields or those experienced during thermal motion could become dominant 
at pressures < 10~^"torr. 
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Figure 7: Photo of the trap with a confined visible graphene flake. The particle is on the left, while scattering 
from the faces of the electrodes is visible on the right. 

7 Measurements 

While the graphene suspensions described above are stable for weeks, for these experiments they are typically 
prepared within 24 hours of their use. After the syringe containing the suspension is placed on the trap column, 
high voltage is applied to the tip, and the gas surrounding it is purged with filtered N2 for several minutes. The 
gate valve is then opened to the trap chamber and liquid flow is initiated out of the electrospray tip by activating 
the syringe pump. Once electrospray emission is initiated, particles arc observable with the video camera 
(typically 0.1-1 scc^^), and pressure in the trap chamber increases. The 500mtorr pressure in the chamber 
used during trapping will decrease, along with collection efliciency, as the pinhole becomes contaminated during 
operation. 

Trapping of particles passing near the trap is facilitated by applying a DC bias (typically ±2 V) to the inner 
electrode that draws the trap minimum towards the electrodes. Once a particle is caught, the gate valve is 
closed; the DC bias can be removed; and the; chamber can be piimpcid to lower pressures with the turbo pump. 
Trapped graphene flakes are highly sensitive to light: at high pressures (when they are trapped) this can be 
seen in "fluttering" (irregular motion) that increases with laser power; at lower pressures (< 10 mtorr) particles 
can rapidly discharge (and leave the trap) when the laser power significantly exceeds 10 mW. These effects are 
consistent with heating of the flakes by absorbed light. Discharging at high light levels occurs at a similar 
rate for both negatively and positively charged trapped flakes, suggesting that thermionic emission is not the 
operative mechanism. Most likely ionic species are present on the flake surface of both polarities and volatilize 
at high temperatures to discharge the flake. 

At the lowest obtainable pressures {p <0.1-1 /xtorr) damping is extremely inefiicient, and particles are fre- 
quently expelled from the trap, most likely from electrical or acoustic noise coupling to the resonant frequency 
of motion of the trapped flake. At higher pressures and at low light levels, particles can remain in the trap for 
weeks. 

During pumping to low pressures, the optical signal coming from scattering from the flake is monitored with 
the photodiode. Strong fluctuations in time appear typically at p < 10 mtorr for graphene flake samples (Fig. 
8a). These fluctuations are largest in magnitude when the sample is viewed with LPL. Such fluctuations do 
not appear when non-graphitic dust particles are trapped. Of the particles that do display fluctuations, greater 
than half of them display virtually 100% signal modulation in LPL at p < O.lmTorr. Only these (~ 100% 
modulating) are chosen for more detailed study. It should be emphasized at the outset that this criterion does 
not mean that the graphene flakes are single layer, since even multilayer flakes are still optically very thin, and 
should display a large anisotropy. It is unlikely, however, that crumpled balls or wads of graphene would display 
large anisotropy, and consequently the measurements discussed below are taken on thin, flat flakes. 

Data for a such a flake at p < 10~^ torr is shown in Fig. 8b. When pumped down while viewing in 
LPL, the fluctuations evident at higher pressures increase to 100% of the photodiode signal. If the sample is 
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Figure 8: (a) Light scattering from a graphene flake at a pressure near the onset of the appearance of fluctuations. 
At higher pressure the scattered Hght signal becomes uniform, (b) Behavior of a flake at low pressure. For the 
circular polarized light measurement, the laser is turned on at i = after being off for an extended period. 
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Table 1: Characteristics of five flake samples. Vtip is the voltage applied during electrospray ionization and 
determines the sign of the charge of the trapped particle. Optical mass is that inferred from Eq. 12 and Fig. 6. 
The rotation resonance measurements were performed on Sample E. 

then illuminated with CPL, the fluctuations almost completely disappear, and the signal approaches a nearly 
constant, equilibrated value. If the illumination subsequently is switched back to LPL, quasiperiodic behavior is 
observed: the characteristic frequency of the fluctuations increases with time until they exceed the measurement 
bandwidth of the photodiode, at which time large amplitude random fluctuations are again observed. 

The likely interpretation of this behavior is that rotation of the flake is induced by CPL, which stabilizes 
its orientation relative to the direction of incident radiation. The quasiperiodic behavior observed when the 
illumination is switched to LPL is a consequence of the gradual slowing of this rotation due to friction with the 
residual gas in the trap chamber. Periodic behavior can occur if there are torques on the flake along an axis 
in the flake plane, which can come from electric fields, as was discussed above. These torques would cause the 
axis of rotation of the flake to precess with angular velocity ~ uj^^/ujg{t). The precession frequency increases as 
the angular velocity of the graphene flake, cjg, slows down, just as the wobbling frequency of a top increases as 
it spins down. 

It has been shown under general circumstances I25J that the axis of symmetry of small, absorbing, oblate 
spheroid will tend to align with the direction of incident radiation when illuminated with CPL. The plane of 
the flake is thus presumably oriented perpendicular to the incident radiation when the particle has reached its 
steady state value when illuminated with CPL. The theory above predicts that the scattered light intensity 
should be maximal in this orientation, since the electric field is always in the plane of the flake. The data in 
Fig. 8b clearly show, however, that peak scattering occurs under LPL illumination and substantially exceeds 
the scattering from the flake after it has been illuminated by CPL for a long period. This discrepancy is 
likely an indication that the flake is not in fact significantly smaller than a wavelength. A flake in the mirror 
configuration, where the angle of incidence is equal to the angle of scattering and where the electric field lies 
in the flake plane, should presumably scatter with peak intensity for a large flake. It is likely that this is the 
orientation of the flake when the signal is maximal under LPL illumination. 

The characteristics of five samples that showed rotating behavior at low pressures are listed in Table 1. q/m 
is determined by observing motional resonance (typically at 1-lOmtorr) and using Eq. 5 with the known trap 

parameters. (^)™pL the peak signal scattered from the sample in LPL, and (^)^'^pL signal from the 

sample after it has been illuminated with CPL for an extended period. Both of these optical measurements are 
made at p < lO-^^torr. Because it is likely to be least sensitive to size effects, (3n)LPL '^^^'^ *° determine 
the "optical" mass using Eq. 12, and q is then determined from the optical mass. 

It is also possible to estimate the mass by observing Brownian motion (also typically at 1-lOmtorr) and 
using Eq. 6. This technique gives results 3-10 times smaller than the optical mass. Systematic errors can easily 
enter this measurement, however: T for the flake under illumination is not known accurately. Also, electronic 
and acoustic noise can contribute to the observed fluctuations and cause the mass to be underestimated. 

Perhaps the most surprising aspect of the data is the extremely long times it takes for the sample to 
equilibrate at low pressures. The rotational damping rates listed in Table 1 are estimated from the changing 
frequency of the quasiperiodic oscillations observed during spin down in LPL, assuming that the torque that 
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is causing precession is constant. The pressure reading on the ionization gauge downstream from the trap was 
between 0.3 and 1.5/xtorr during these measurements. The resuhant values for p/^^ are large and suggest that 
the samples all contain many layers if the simple kinetic picture of damping presented above is correct and if 
lu> ~ 7t)- The pressure measurement in the experiment is certainly inadequate, but to bring the data into line 
with theory for a monolayer flake, the pressure in the chamber would need to be lower than the gauge value, 
which would be surprising. (The cold cathode gauge reading is closely matched to the capacitance diaphragm 
gauge, very near to the trap center, when they are both in range.) 

8 Observation of Rotational Resonance at MHz Frequencies 

The data presented above provides strong indirect evidence that flake spinning is induced by exposure to CPL. 
The long measured values of - combined with the estimates of optical torque for graphene - suggest that the 
rotation frequency is in the MHz range in the experiments. Direct optical observation of this spinning using 
fast optical detectors is a challenge due to the small (femtowatt) signals scattering off the flakes. An alternative 
is to excite the flakes at high frequencies and seek to observe changes in the optical signal measured at low 
frequencies. In particular, a sample can be exposed to high frequency electric fields by applying a voltage to 
the inner electrode of the trap. At MHz frequencies, these fields have a negligible effect on the position of the 
trapped fiake, but do produce torques on the sample as was discussed above. 

To perform this experiment a sample was illuminated with CPL at p < 10 /xtorr. An 8 V amplitude AC 
signal was applied to the inner electrode, corresponding to a peak E at the trap center of 1250 V-m^^. The 
frequency was ramped logarithmically while the optical scattering signal was monitored (Fig. 9). The data 
clearly show sharp jumps in the photodiode signal, occurring both when the frequency is scanned upwards and 
downwards. The jumps are not uniform in size, however. Note also that the direction of the jump (towards 
brightening or dimming) is not always the same for a given direction of the frequency sweep. 

The rotation resonance jumps show clear dependence on the pressure in the trap: the magnitude of the 
jumps decreases markedly at low pressure, and the frequency at which the jumps occur is roughly inversely 
proportional to the pressure (Fig. 10). The scatter in the frequency associated with the jumps is increasing 
at high pressures, probably an indication that random or thermal fiuctuations in the rotation and orientation 
of the fiake are becoming significant. Subsequent to these resonance measurements, the bias applied to the 
inner electrode was switched off, and was determined using observation of the spin down behavior in linear 
polarized light (Fig. 11). 

The fact that signal intensity changes abruptly strongly suggests that the axis of orientation of the fiake (or 
its average orientation) is shifted by a resonant interaction with the applied torque. While it is tempting to 
treat this as a problem of rigid body dynamics, a better approximation for graphene is to examine the dynamics 
of a thin membrane with no rigidity whatsoever under tension purely from centrifugal forces. This problem was 
treated by Lamb and Southwell [4^ for the case of a circular disk. The lowest frequency transverse excitations 
of the disk spinning with angular velocity uj have a single nodal diameter that rotates at ±0; in the frame moving 
with the spinning disk. In a nonrotating frame, these excitations correspond to a DC tilt of the axis of rotation 
and fast precession at 2a; of the rotation axis of the disk. This latter motion can presumably be excited by an 
external torque applied at angular frequency 2cj. 

For rigid 2D plates of arbitrary shape (where + 1\\2 = I-l) i small deviations of the instantaneous axis of 
rotation from the principal axis perpendicular to the plate precess at angular velocity [45] : 



in the co-rotating frame of reference, just as was the case for the spinning circular membrane. It is likely 
then that under very general conditions, rotating 2D objects are susceptible to transverse excitation at angular 
frequency 2a;. 

Equating the observed resonance frequency in the data with twice the rotation frequency of the flake, ujg/n, 
allows for the estimation of the flake dimensions: "f^ may be estimated from the data in Fig. 11. Combining these 
values (both measured at 1.5/itorr) allows for the determination of uj± = a;g7^ = 1.7 x lO^sec"^. Eq. 16 can 
then be used to evaluate the size of the fiake: a =0.4 /im. Furthermore, p/^^i can be used to estimate n=17, if the 
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Figure 9: Observation of rotational resonance of a graphene flake. For all of the data, the sample is exposed to 
4niW circular polarized light during the measurements. 
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Figure 10: Effect of chamber pressure on the frequency at which rotation resonance jumps occur. 



simple kinetic theory of damping is vahd. A flake with these dimensions would have m = Tra'^np2u = 6 x 10^^^ kg, 
a factor of four lower than the optical mass listed in Table 1. 

Finally, using these estimates for the particle size and the result of Eq. 20, the precession rate of the 
axis of the flake induced by torques from the resonant electric field can be estimated: uj^^/cog '-^0.3sec~^(at 
p=1.5 /xtorr). Resonant reorientation (on the slow experimental time scales) is thus a reasonable explanation for 
the observed modulation of light scattering from the rotating graphene flake at specific excitation frequencies. 
A detailed understanding of the reorientation dynamics observed in Fig. 9, however, will require additional 
exp crime nt at ion . 

9 Optimizing the Measurements 

The measurements presented above are obviously preliminary, and many aspects of the experimental design 
need to be improved. More uniform samples and better techniques to estimate their dimensions are desirable. 
A more accurate model of light scattering from micron-scale graphene flakes will also be required. Below is a 
list of several improvements to the experiment. 

9.1 Trap Design 

The trap design presented in Fig. 1 suffers from the disadvantage that there is a large cubic term in the confining 
pseudopotential near the trap minimum. When excited to large amplitudes the system exhibits hysteresis, 
making it difficult to measure the frequency of resonance accurately. Traps with a symmetric design can make 
much more precise determinations of q/m ^1, and (by observing single electron discharging or charging events) 
of q and m separately |22| . An advantage of the current design, however, is that a trapped particle can in 
principle be transferred between two traps whose apexes are brought close together. Sample transfer may be 
useful, for example, in a load lock to a UHV or cryogenic environment. 

9.2 Particle Injection and Preparation 

While the techniques presented above for creating graphene suspensions and injecting them into the trap have 
yielded encouraging preliminary results, the data suggest that that flakes trapped so far have been multi- 
(or even many) layer. It is possible that the suspension contains mostly multilayer flakes. However, it is 
also likely that there is substantial selection bias in the experiment towards trapping large particles: given a 
^30 mg-L~^suspension of graphene injected into the trap chamber at 1 /iL-min~^, there are 10^ — 10^ particles 
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Figure 11: Spin down behavior measured subsequent to observation of rotational resonance of the flake observed 
with 4niW of hnear polarized light. Previous to i = 0, the flake had been exposed to 4mW of circular 
polarized light for a long period. To facilitate the observation of the quasiperiodic oscillations, the time axis is 
exponentially expanded in the bottom plot. 
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injected every second, while typically ^0.1-1 sec~^are observed to enter the trap vicinity. If the trap design 
favors trapping heavy particles (for example if the chamber pressure is so high that the large damping of light 
particles prevents effective trapping), then measurements of trapped particles will not reflect the distribution 
of sizes of particles in the original suspension. Further measurements on both the distribution of particle sizes 
in the suspension and of biases present in the trapping design are necessary. 

In the current experiment q/m is not a readily adjustable parameter (aside from its sign). Further charging 
of trapped particles would be desirable, for example, to offset the discharging of the particles that was observed 
at high light intensities or to increase flake charge beyond what is possible using the electrospray ionization 
method. Additional charging could be provided, either by using ultraviolet light [H] or by an electron beam 
[22]- 

Finally, while suspensions of graphene created using ultrasonication yield irregularly shaped flakes with a 
distribution of sizes, substantial improvement are possible: with a starting material of long cylinders etched 
using deep reactive ion etching (RIE) [35] from highly oriented pyrolytic graphite (HOPG), for example, the 
lateral dimensions of the flake could be made uniform. Ultra-centrifugation techniques |47) may lead to the 
production of suspensions containing flakes of only a specified number of layers. 

9.3 Optics 

The optical setup presented above will clearly benefit from a more rapid and versatile method for controlling light 
polarization. The measurements of rotational resonance were performed under continuous CPL illumination 
for simplicity. A better experiment would illuminate with CPL to confer angular momentum on the flake and 
measure with LPL in order to maximize sensitivity to flake orientation. Another improvement would be to view 
scattering at an angle perpendicular to the incoming radiation (and thus in the plane of the spinning disk) to 
maximize sensitivity of scattering to small changes in the orientation of the flake. 



10 Applications for Levitated Graphene 

10.1 Mechanical and Materials Properties Measurements 

Graphene has been deemed the world's strongest material because of its large Young's modulus and proven 
ability to withstand tensile strain in excess of 10% [11]. Spinning graphene up to high rotational velocities 
should provide an important measurement tool for graphene in an environment where uniform and controllable 
tensile stresses can be applied. The tensile stress from centrifugal force, fc, at the center of a circular graphene 
single layer rotating about an axis perpendicular to its plane is [H]: 

/c= ^(3 + i^)aWp2D, (23) 

where f is the Poisson ratio (around ~0.17 [IT] for graphene). If the yield strength is roughly at 10% strain and 
the Young's modulus is i?2D=340N-m~^[TT], then the maximum possible value of aoj ^ 10^m-sec~^, and the 
maximum rotation frequency of a graphene flake with a—1 fim is ~1.7GHz. While this rotation rate is much 
greater than was achieved in the experiments presented above, higher rotational velocities should be attainable, 
either by performing experiments at lower pressures or at higher laser powers. 

In addition to centrifugal force, charged graphene also experiences electrostatic tension, /g, which may be 
estimated from the capacitance (C) of a circular disc [43] : 

/ (24) 
■'^ d{iTa^)\2CJ 27ra da V 16eoa y 327reoa3 ' ^ ' 

For the experiments discussed above, q ^ ±2000 e and a ~0.5/im so ^ 10^'^N-m^^. While this number 
suggests that flake charge can be substantially increased, it is likely that electronic or ionic field emission at 
the edges of the flake will determine the maximum charge it can hold, rather than the intrinsic strength of the 
graphene. 

The measurement of rotational resonance, in addition to providing direct information about the flake's 
rotation frequency, may provide a means to measure dissipation of vibrational excitations of the spinning disk 
when coupling to the external environment is extremely small and controllable. Rotational resonance may 
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have sensitivity exceeding mechanical resonance of graphene attached to substrates for determination of the 
mechanical properties of the material, especially in the regime of large tensile deformations. 

The thermal isolation of trapped graphene, combined with its strong optical absorption, means that it should 
be relatively easy to measure graphene's material and chemical properties at extremely high temperatures 
without any interaction from a substrate. Very little is known about the melting of graphene or how defects 
will behave at high temperatures. Temperature dependent adsorption and desorption of various atomic species 
|48j introduced into the trap chamber can presumably be measured with high accuracy, either by measurement 
of the flake mass [55] or by rotational resonance. 

10.2 Crystal Growth and Modification of Trapped Graphene 

A significant limitation of the approach to graphene trapping presented above is that flake size is limited: 
it is improbable that flakes larger than ^10 fim can be injected into the trap using liquid suspensions and 
the electrospray technique. Trapped flakes may possibly be modified in situ, however. For example, high 
temperature anneals of rotating flakes might alter the shape of the trapped flake by promoting migration and 
smoothing rough edges. It is possible that the trapped environment may also have applications for graphene 
crystal growth: the general arguments against 2D crystallization HHj are unlikely to be valid in the presence of 
electrostatic or centrifugal tensions that would tend to keep the structure planar. The low mass of graphene 
monolayers means that even "wafer-scale" graphene can be held in traps without difficulty. The challenge is to 
find appropriate conditions where a small trapped graphene crystal will expand and maintain crystallinity when 
exposed to carbon sources, such as those used for C doping in molecular beam epitaxy [Sn][5T]. One possibility 
is that vacancies are first injected into the bulk from the edge by a brief exposure to a high temperature. 
Subsequently, at a lower temperature, C from an impinging molecular beam is incorporated into the crystal at 
the positions of the vacancies. Because of the short time constants associated with cooling graphene at high 
temperatures (Fig. 5), such multi-temperature growth cycles could proceed very rapidly. 

The fact that electrostatic or centrifugal tensioning of graphene flakes can be used to increase the lattice 
constant by up to 10% may also facilitate novel heteroepitaxial materials based on a graphene substrate with an 
adjustable lattice. Boron nitride, with a lattice constant about 2% greater than graphene, may grow conformally 
on an appropriately tensioned substrate. Possible C-BN multilayers grown in this way may modify the band 
structure of the carriers in the graphene |52j or improve their mobility 7 . 

Is it likely that for many applications a sample prepared or modified in a trap will ultimately need to 
be positioned on a substrate. It is possible that electromagnetic focusing techniques (like those used in mass 
spectroscopy and electron microscopy) could be used with trapped flakes to place them accurately on substrates, 
although maintaining proper orientation of the flake during its deposition onto a surface could be challenging. 

10.3 New Physics 

The strong optical absorption that makes observing and spinning graphene easy will unfortunately make if very 
challenging to cool samples and observe new physics, either of the mechanical system or of the confined charge 
carriers. The optical measurement scheme presented above would need to be dramatically improved to enable 
low temperature measurements. Nonetheless it is worth mentioning some experimental possibilities worthy of 
investigation: buffer gas cooling of the trapped flake using •^He may enable cooling to ^0.3 K. Measurement of 
the 2D charge carriers could be accomplished either by their effect on the rotational resonance behavior or by 
the braking effect noted above (Eq. 21). Graphene samples annealed at high temperature to remove adsorbates 
and measured in high vacuum may have extremely high carrier mobilities [3] . Traps could be placed in magnetic 
flelds to facilitate measurements in the quantized Hall regime, but the charge density on the flakes will not be 
uniform [33] and will be sensitive to the orientation and the rotational velocity of the flake. Torques associated 
with orbital diamagnetism [53] of the 2D charge carriers should be detectable by their effect on flake orientation 
|54j . It is also possible that under appropriate conditions resonant coupling between electron or nuclear spins 
situated in the flake and the rotation or orientation of the flake could be observable. 

Finally, it may be possible to cool levitated graphene using unorthodox techniques: for example, adiabatic 
detensioning of the vibrational modes of a spinning membrane may lead to cooling as the rotational frequency is 
slowed. Low temperatures could conceivably be reached either by coupling to a laser-cooled atomic system |55j 
[56] or by direct optical cooling using cavity optomechanics techniques [57] [58] [59]. Perhaps the ultimate goal 
is to reach kT/h < Wg, a regime where thermal excitations of the rotating membrane are suppressed, analogous 
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to cooling to the ground state of a mechanical system |60j . For a graphene flake rotation frequency of order 
1 GHz, this goal can be reached at T ~ 50 mK. Trapped graphene would then be a "spinning qubit" for studies 
in quantum information science. 
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